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Abstract. We consider a cell population described by an age-structured par- 
tial differential equation with time periodic coefficients. We assume that divi- 
sion only occurs after a minimal age (majority) and within certain time inter- 
vals. We study the asymptotic behavior of the dominant Floquet eigenvalue, 
or Perron-Frobenius eigenvalue, representing the growth rate, as a function of 
the majority age, when the division rate tends to infinity (divisions become 
instantaneous). We show that the dominant Floquet eigenvalue converges to 
a staircase function with an infinite number of steps, determined by a discrete 
dynamical system. As an intermediate result, we give a structural condition 
which guarantees that the dominant Floquet eigenvalue is a nondecreasing 
function of the division rate. We also give a counter example showing that the 
latter monotonicity property does not hold in general. 



1. Introduction 

1.1. Age structured model of the cell division cycle. Among all the popula- 
tion models, age-structured population models are perhaps the simplest and 

oldest ones taking into account the variability between individuals. Age struc- 
tured partial differential equations have been introduced in epidemiology through 
the famous Von Foerster - Mac Kendrick model [] 2] (where age stands for age in 
the disease). In its simplest version, it describes the dynamics of a population rep- 
resented by a density n(t, x) of individuals of age between x and x + dx through 
age-dependent birth and death rate: 
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( dn dn T , . . . 
-^ + -^ + d(x)n(t,x) =0, 



CD 



n{t, 0) 



B(x)n(t, x)dx, 



, n(0, x) = n (x) given. 



A special case of this model can be used to represent cell division, taking 
d(x) = K(x) and B(x) = 2K(x). More complex models, in terms of structure, 
representing the cell division cycle (linear or not, size or cyclin structured) have 
been developped and studied in [16, 10, 21]. In this paper, we study an extension 
of the age-structured division equation, in which the coefficients are periodic func- 
tions of time. There are several motivations to introduce time dependence in (1). 
The most classical one is the representation of seasonality (see [2] for instance). 
Concerning cell division, our main motivation comes from cancer therapy such as 
resonance therapy [8] or chronotherapy [13]. The latter is based on the use of circa- 
dian rhythms and motivates especially the introduction of time periodic coefficients 
(see [7, 11] for instance). We also refer to [9, 14, 19] for work on renewal models 
with periodic coefficients. 

In our case, the periodic division equation reads 

dn dn 
dt dx 



(2) 



n(t,0) -- 
n(0, x) 



K(t, x)n(t, x) = 0, 
K(t, x)n(t, x)dx, 



/o 

n°(x) 



given. 



where K satisfies K(t + T, •) = K(t, •) for some T > 0. 

For equations (1),(2), the population grows exponentially with time (see [18] for 
instance) and the growth rate is measured by the so called Malthus parameter A 
(n(t,x) ~ pe xt N(t,x), N being given by (3)). Particularly, in (2), the Malthus 
parameter (F for Floquet) can be generally described by means of a Floquct 
eigenvalue problem. 



(3) 



ON 


dN 


1H + 


dx 


n(t,0) 


= 2 


N(t + 


T,x 



(X F +K(t,x))N(t,x) = 0, 

X) 

K(t,x)N(t,x)dx, 



together with the dual eigcnproblcm 
(4) I ~di ~ dx 

{ 4>(t + T,x) = 4>(t,x), 4>>o 



+ {X F + K(t, x))(j>{t, x) = 2K(t, x)</)(t, 0), 



For uniqueness issues, these equations are usually completed with the following 
normalization: 



(5) 



o Jo 



Ndxdt 



o Jo 



N(t,x)(j)(t,x)dxdt = 1. 
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Actually, the integral J °° N(t,x)<fr(t,x)dx does not depend on t. Various proper- 
ties of the Floquet eigenvalue have been studied in [2, 3, 4, 5, 6] with in mind the 
comparison with time independent systems (in which K (t, x) is replaced by suit- 
able time averages). Here, we are interested in properties of the eigenvalue (and 
more generally of the dynamics) that are intrinsically related to the fact that the 
coefficients are time-dependent. 

We consider in Section 3 a division rate of the form 

(6) K(t,x) = KtP(t)B(x)l [ai+oo[ (x) . 

Here, k > is a scaling parameter and the function tp(t) ^ expresses the de- 
pendence of the division rate as a function of time. We denote by 11 the indicator 
function of a set. The term S(x)lr Oj+00 [(x) indicates that the cells only divide 
when x ^ a, i.e., when they are older than a "majority age" a. The function B(x) 
expresses a general modulation of the division rate as a function of the age. 
The function B will be required to be bounded, positive, of infinite integral, 

(7) B E L°°(R + ), B > 0, \ B = +oo 

This is needed for the existence theory of the Floquet eigenvalue (see the appendix) . 

1.2. Main results. Our goal is to derive the asymptotics of the growth rate when 
the division becomes instantaneous, meaning that the scaling parameter k in (6) 
tends to oo. The growth rate is represented by the Floquet eigenvalue, Xf, which is 
the unique scalar A solution of (3), (4). Hence, we shall investigate the asymptotic 
behavior of Xf as a function of k. 

To this end, we first study the monotonicity of the Floquet eigenvalue with 
respect to the division coefficient K(t, x). Since every division replaces one old cell 
by two new ones, it is natural to think that the growth rate is a nondecreasing 
function of K{t,x). Our first result is an explicit counter example, showing that 
this apparently intuitive property is not valid due to the time and age dependence 
of the division rate. 

Proposition 1. There exists configurations such that 

V(t,x), < K 1 ^) < K 2 (t,x), 
and the corresponding Floquet eigenvalues satisfy 

Xp > Xp. 

However, the following theorem identifies a subclass of transition rates for which 
the Floquet eigenvalue is indeed a monotone function of the transition rate. 

Theorem 1.1 (Comparison principle). If K 1 satisfies 

v n> / K (s,s — v)ds is nondecreasing for any t, 

J V 

then for all K 2 > 0, 

(V(*,ar) K 2 (t,x) ^ K%x)) Xp ^ Xp, 

(V(t,x) K 2 (t,x) < K^t^x)) => X F ^ Xp. 
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The proof relies on a disaggregation idea of independent interest. We first show 
that the density n(t,x) can be written as ^2 i>0 2 l m(t,x), where every rii(t,x) can 
be interpreted as the density of an individual of the zth generation, and the densities 
(^i)i^o satisfy an infinite triangular system of transport equations coupled by inte- 
gral terms. Then, for all j ^ 0, we show that the term Sj(t) = J2i^j Iq° n i(*> x)dx 
is a nondecreasing function of K S {K 1 ,!^ 2 }, and derive from this fact the com- 
parison of Floquet eigenvalues. We note that the latter fact has a probabilistic 
interpretation: it means precisely that the generation of a single individual is a 
nondecreasing function of K G {K 1 , K 2 } with respect to the stochastic (majoriza- 
tion) order. 

When inf ip > 0, the intuition predicts that the growth rate goes to the limit 
log 2/ a as k tends to infinity, because every cell will divide shortly after reaching 
age a when the parameter k is large. Our next result confirms that this is the case. 

Proposition 2. Suppose that tjj is a bounded T -periodic junction with w£ip > ; 
and that B is a positive bounded function satisfying (7). 
Then, for all k > 0, 

(8) <A F (*)< 

a 

and 

(9) lim \f{k) = 

k— >oo a 

When ip vanishes, meaning that the division of cells is blocked at certain times 
of the day, we shall see that complex synchronization phenomena appear. We shall 
assume that ip has a square wave shape, i.e., that ip is a T-periodic function such 
that, for some < r < T, 



(io) m = 



1 for t e [0, t) 
for t e [r, T) 



The following is the main result of this paper. We denote by \x] the smallest 
integer which is not less than a real number x. 

Theorem 1.2 (Discrete limit). Suppose that the division rate is given by K(t, x) = 
K'0(t)]l[ aiOO )(x) ; where ip is given by (10). Denote W T = [r, T] + NT and let N a be 
the integer defined by 

N a a = min W T CiNa . 
Then, the limit of the Floquet eigenvalue A K (a) as k — > +oo is given by 

(11) A°»= ^ log2. 

\N a a/T\T 

Moreover, we have the estimate 

(12) A°°(a, r) ^ A K (a, r) ^ A°°(a) ^1 - e~ Kr ^ T ^ , 
where the function r(a, r) (the convergence rate) satisfies 

(13) r(o,r)> min h^, T \ a r = sup{a', \°°(a') = A°°(a)}, 
and r is positive outside a closed countable set. 
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The function a h-> A°°(a) is a staircase with a countable number of steps (see 
Figure 1 below and Figure 8 and Figure 9 in Section 3.5). The quantity a r — a which 
controls the rate represents the distance to the right end of the step containing a 
(on Figure 1 for instance, we can observe that the speed of convergence depends on 
this distance). 



— K=10 

— K=30 
K=50 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

minimal age lor division (a| 

FIGURE 1. Convergence for t = 2/3 and T = 1. 

There is a special case in which the limit of the Floquet eigenvalue can be rewrit- 
ten in terms of integer parts. If t/T < 1/2, and < a < T, then, we have 

A°°(o)= r-llog2 . 
a 

We shall derive intuitively Formula (11) in Section 3.2.1, by introducing a discrete 
dynamical system, equipped with a multiplicative functional, the asymptotic geo- 
metric mean of which yields A°°(a). The existence of such a simple asymptotic 
formula reflects the fact that the Perron-Frobcnius operator describing the evolu- 
tion of the population degenerates as k — > oo. The rate r(a) appears naturally in 
the proof of convergence. This proof relies on Theorem 1.1 (the comparison prin- 
ciple), on the generational model already used in the proof of this theorem, and on 
explicit bounds of the density, for a special initial condition leading to a trajectory 
"close" to that of the discrete dynamical system. 

The convergence rate in Theorem 1.2 is reminiscent of large-deviation or Laplace 
asymptotics. However, this theorem differs in its essence of already known asymp- 
totics results in Perron-Frobenius theory with such large deviations flavor, includ- 
ing zero-temperature asymptotics of Ising type models or Frcnkcl-Kantorova mod- 
els [1, 20]. In these problems, the limit of the Perron eigenvalue, in a log-scale, 
turns out to be the ergodic constant of a controlled deterministic dynamical sys- 
tem, whereas in the present case, the dynamical system which determines the limit 
involves no control. 

When ip vanishes, we may consider t/j £ (t) := ip(t) + e, with e > 0, and define 
the growth rate A(e, n) with an obvious notation. A comparison of the previous 
theorems shows that the limits in k and e do not commute 

= lim lim X(e,K) + lim limA(e,/c) = A°°(a) . 
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2. MONOTONICITY WITH RESPECT TO THE DIVISION RATE 

Wc first focus on the question of the monotonicity of the Floquct eigenvalue with 
respect to the division rate. 

2.1. A counter example in which the growth is not a monotone function 
of the division rate. Wc next show that increasing the division rate may decrease 
the growth rate of the population. The counter example is based on the following 
exactly solvable model. We assume that the period is T — 1. and fix a G (0, 1). It 
will be convenient to write 

h = [0, a) + Z, h = [a, 1) + Z , 

and 

Xj {t) := llj. (t), j = 1, 2, so that X i + Xi = 1 • 
We shall assume that the division rate is of the form 

K(t,x)= X i{t-x)K 1 (t)+X2(t-x)K 2 (t) . 
Integrating along a characteristic, we get 

n(t, x) = n{t — x, 0) exp ( — / K{t — x + y, y)dy) 

Jo 

and so 

n(t, x) = n(t — x, 0) exp ( — / Kj (t — x + y)dy) , if t — x € Ij . 

Jo 

In other words, when an individual is born at time t S Ij, then its next division 
will occur with a rate Kj(t). This model has an intuitive interpretation. We may 
assume for instance that I\ represents the days, whereas I2 represents the nights. 
If being born during the night or during the day influences the fertility (division 
rate) of the individuals, we arrive at this model. We may then consider that there 
are actually two subpopulations rij(t, x), j = 1,2, corresponding respectively to 
individuals born during the day or during the night. Wc have 

nj(t,x) =n(t,x)xj(t- x), 3 = 1,2 

They satisfy the equations: 

d t rij + d x nj + Kj(t)rij = 0, j = 1,2 

/>oo 

n j (t,0)=2 Xj (t) (K 1 (t)n 1 (t,x)+K 2 (t)n 2 {t,x))dx j = 1,2 . 



Denoting by 

= / rij(t,x)dx 



the total population of type j at time t, we get that the vector P(t) := (P\{t) , P 2 {t)) T 
satisfies the ODE: 

|p(t) = M{t)P{t) 

where 

' (2xi (i) - l)K x {t) 2 X i(t)K 2 (t) 
2 X 2(t)K 2 (t) (2 X2 (t) - l)K 2 {t) 



M(t) 
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We now assume that the functions K\(t) and K 2 (t) are constant during the day 
and during the night, meaning that 

Kj.it) = oixi(t) + 6ixa(t), K 2 (t) = a 2X i{t) + 62X2 (t), 

where the coefficients o»,6j are constant. With this form of the coefficients, we 
obtain 

M(t)= X i(t)M a + xn(t)M b 

where 

M " = ( ° - ) ' * = ( * * ) ' 

An important point here, still with the above interpretation is the following: a 2 and 
b\ could be considered in some sense as transition coefficients between populations 
Tii and n 2 . 

Since the linear dynamics M(t) switches between the constant dynamics M a and 
M b , which are successively exercised during the intervals [0,a) and [a, 1), we get 
that 

P(l) = exp ((1 - a)M b ) exp (aM a )P(0) . 

Denoting by p the Perron root (spectral radius) of a nonnegative matrix, we arrive 
at the following expression for the Floquet eigenvalue 

(14) X F = log p^exp((l- a)M b ) exp (oM„)j . 

Let us now consider the situation in which a 2 = b 2 = 0, meaning that an individual 
born during the night is sterile (division will never occur) . Then, the second columns 
of the matrices M a and M b vanish, and so 

/, > n / n / exp(aai - (1 - a) 61) 0\ 
exp ((1 - a)M b ) exp (aM a ) = I ^ v ^ 1 J 

(the value of the off diagonal entry, denoted by *, is irrelevant). Since the spectral 
radius of a nonnegative triangular matrix is the maximum of its of its diagonal 
entries, we get 

\f = max (aai — (1 — a)6i, 0) . 
Hence, is an increasing function of the division rate ai, but a decreasing function 
of the division rate b\ . This establishes Proposition 1 . □ 
It should be noted that the preceding counter-example subsists when a 2 and 
62 are sufficiently small, because the spectral radius is a continuous function of 
the parameters. This is confirmed by Figure 2, in which the Floquet eigenvalue is 
plotted as a function of b± and b 2 for a\ = 10 and a 2 = 0.1. 

Remark 1. This counter example can be understood intuitively in Figure 3. The 
values of the division rate K(t, x) are represented (Left, the values are constant on 
each cell). Consider now a chain of descendants of the same individual (i.e., this 
individual, one of his children, one of the children of this children, etc.). This can 
be represented by a Markov process X t with state space in K + (the age), which 
jumps from age x to age with rate K(t,x). Each time the path hits the axis 
x = corresponds to a division. If the path hits x = during the night, since 
b 2 = 0, the corresponding individual is sterile (no division will ever occur, which is 
reflected that the path is now an unending line of slope 1 in the (i, x) plane). The 
production of sterile individuals occurs when the process jumps from a cell labeled 
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b1 



Figure 2. The counter example. The Floquet eigenvalue as a 
function of the division rates b\ and 62, when a\ — 10 and 02 =0.1. 
When 62 is small, increasing b\ (and so, increasing the division rate 
K(t,x)) decreases the growth rate. 

"61" to the x axis, and this occurs precisely with rate b\, which explains intuitively 
why the growth rate, A^, is a decreasing function of b\. 




a I t t 



Figure 3. The counter example explained. The values of K(t,x) 
(left). The night is shown in grey. A typical path (right, broken 
line). The child born during the night is sterile (final segment of 
the path), and the division coefficient b\ determines the creation 
of this child. 

Remark 2. The system does not enter in the framework considered in Appendix 
for the existence of the Floquet eigenvalue. Indeed, the preceding study merely 
determines the growth rate of the "aggregated" population P(t), which is deter- 
mined as the Perron eigenvalue of a matrix. However, the growth of n(t, x) is 
readily derived from the one of P(t). To see this, it is convenient to introduce the 
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adjoint positive Floquet eigenvector of M(t), whose existence and uniqueness are 
guaranteed by Perron- Frobenius theory (as soon as the coefficients a\, 0,2, &i, 62 are 
positive): 



dt 

We have 



^ ( fa fa ) + A ( fa fa ) + ( fa fa) M(t) = 0. 



f t (fat),P(t)) = X(fat),P(t)) 



which can be read as, denoting fat, x) = fa(t)xi(t — x) + fa(t)\2(t — x), 



— I n(t,x)fat,x)dx = A / n(t,x)fat,x)dx, 
Jo 



dt 



It is straightforward to check that fa and fa are positive (still assuming that the 
coefficients 01,02,61,62 are positive). Hence, there exists M ^ m > such that 
m ^ fat, x) ^ M, and therefore, 

me xt / n°dx < / n(t,x)dx ^ Me xt / n°dx. 



This shows that the Floquet eigenvalue A computed from the aggregated dynamical 
system does determine the behavior of the solutions. 

2.2. A sufficient condition for monotonicity. The preceding counter example 
is due to a "birth day penalty" (the date of birth of the individual influences 
critically its division rate) . We shall avoid such a pathological behavior by making 
the following assumption: 

(15) (t, v) 1— > / K(s, s — v)ds is nonincreasing in v. 

J V 

This condition is fulfilled when the division rate K(t, x) ^ depends only on one of 
the two variables t and x. It is also fulfilled if for any i, K(t, •) is a nondecrcasing 
function (and particularly, in the case of separated variables K(t,x) = ^i{i)B{x) 
when B is nondecreasing). 

This is illustrated in Figure 4. The condition requires the integral of K over 
the characteristic shown on the figure to be nonincreasing in v. Note that this has 
a probabilist interpretation in terms of the Markov process X t introduced in the 
preceding section. Indeed, 

exp ( — / K(s, s — v)ds) 

J V 

represents the survival probability at time t of an individual born at time v, i.e., the 
probability, conditional to X v = 0, that the process makes no jump to point before 
time t. The condition requires this probability to be a nondecreasing function of v, 
in other words, that by delaying the birth, the survival probability at a given time 
cannot decrease. 

Theorem 2.1. Provided K 1 or K 2 satisfies (15) and K 1 ^ K 2 , then, for any 
n° G L 1 (M + ), if we denote by n l the solution of (2) with K = K l , we have, for any 
t ^ 0, 

/>oo />oo 

n 1 (t,x)dx^ / n 2 (t,x)dx. 
Jo 
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t rt 

K(s,s — v)ds^ / K(s,s — u)ds 



(16) 



Figure 4. The sufficient monotonicity condition: the survival 
probability at time t (the exponential of the opposite of the in- 
tegral) increases with the birth time v. 

To establish this result, we introduce the following system of PDEs: 
' dtUi + d x rii + K(t, x)rii = 0, i > 0, 
no(*,0) =0 

/>oo 

rii(t,0) — / K(t,x)rii-i(t,x)dx, i^O 
Jo 

, no(t = 0, x) = n (x), rii(t = 0, x) = 0, for i ^ 1. 



We easily check that n = J^^o^ 71 * ^ s the solution to (2). Intuitively, the 
term 2 l rii represents the number of individuals of the "generation" i, and so, we 
shall refer to this model as the "generational model" in the sequel. We denote by 
(nj)i^o, (resp. (^f)i^o) the solution to (16) with K = K 1 (resp. K = K 2 ). We also 
introduce: 



/•OO 

S j(t) "•= X) / rn{t,x)dx, for j ^ 
Jo 



A short computation leads to 



/>OC 

/ n{t,x)dx = S Q (t) +V2 J - 1 S J (() 

7 ° 3 >1 



Lemma 2.2. VFe /aave: 

poo 

(17) S (*) = 5*o(0) = / n°(x)dx , 

Jo 

(18) Si(t)=J^ n°(x)(l-exp(- K(s, x + s)ds) S j dx , 
and /or j ^ 1, 

(19) Sj{t) = ^ n j - 1 (s,0)(l-exp(-J K(s + y,y)dy)^ds 
Proof. It is straightforward to show that for every j ' ^ 0, 



-5 J (t)=n i (t,0) , 
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and so 

(20) Sj(t)= [ nj(s,0)ds + Sj(0) . 

Jo 

This leads immediately to the first statement, since for all s > 0, no(s, 0) = 0. 
To establish the next statements, we shall use the characteristics: 

no(t, x) = 0, if x ^ t, 

(21) n (t, x) = n°(x — t) exp ( — / K(s, x — t + s)ds), if x ^ t , 

Jo 

and for i ^ 1, 

rii(t, x) = 0, if £ ^ .t, 

(22) Uj(i, a;) = 7ii(i — x, 0)exp ( — / K(t — x + y,y)dy) , if x ^ t . 

Jo 

The lack of symmetry between the cases i ^ 1 and i = is due to the boundary 
conditions: no(0, x) = n°(x) is known, whereas for i ^ 0, n.i(i, 0) can be inductively 
assumed to be known, using n\-\. 

Since n,i(0,x) = 0, we have <Si(0) = 0. Using (20), we get 

/>t />t />oo 

1S1 (t) = / ni(s,0)(is= / / K(s,x)no{s 7 x)dxds 
Js=o Js=o Jx=o 

ft />00 />S 

= / / K(s,x)n° (x — s) exp ( — / K(u,x — s + u)du)dxds 

J s=0 J x=s Jo 

' K (s, y + s)n°(y) exp ( — / K(u,y + u)du)dyds , 

l s =0Jy=0 Jo 

and, interchanging the order of integration, 

Si(t) = / n°(y) / K(s, y + s) exp ( - / K(u,y + u)du)dsdy 
Jy=o J s=o Jo 

f n°(y)^l-exp(-^ tf(s, y + s)ds)^ dy, 

which is the second statement. 

Assume now that j ^ 2. Arguing as above, but using this time (22) instead 
of (21), we get 

r-t poo 

= / / K(s,x)rij-i(s,x)dxds 

J s=0 Jx=0 



s=0 Jx=0 



K(s, x)rij-i (s — x, 0) exp ( — / K(s — x + y,y)dy)dxds 

Jo 

K(s,x)n,j_i(s — x,0) exp ( — / K(s — x + y,y)dy^jdsdx 



x—0 J s—x 

t pt—X pX 

/ K (u + x, x)rij-i(u, 0) exp ( — / K{u + y,y)dy)dudx 

x=0 Ju=o Jo 

t pt — U px 

rij-i(u,0) / K(u + x, x) exp ( — / K(u + y,y)dy)dxdu, 

=o Jk=o Jo 

which leads to the third statement. □ 
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Theorem 2.1 will be obtained as an immediate consequence of the following 
lemma. 

Lemma 2.3. Suppose A 1 satisfies (15), then 

A 2 ^ A 1 =► S]{t)>S}(t), Vt,j 

K 2 < K 1 =► S](t)^Sj{t), Vt,j . 

Proof. The proof is performed by induction on j. We shall only consider the case 
in which K 2 A' 1 , the other case being similar. 

For j = 0,1, the conclusion follows readily from the two first statements in 
Lemma 2.2. Assume now that j ^ 2. We set 

i*t — V 

Q k (t,v):=l-exp(- K k (v + y,y)dy) 



Using (19), we get 



S}(t)-S](t)= / n)_ x {v, 0)Q 1 (t, v)dv - / n 2 j _ 1 (v,0)Q 2 (t,v)dv, 



{n)_ x {v, 0) - n 2 (v, 0))Q\t, v)dv + / n 2 {v, 0)(Q 1 (t, v) - Q 2 (t, v))dv 



n 2 _ 1 {v,0){Q 1 {t 7 v) -Q 2 (t,v))dv, 

= -J (Sj_ 1 (v)-S 2 _ 1 (v))^Q 1 (t,v)+ J n 2 _ 1 (v,0)(Q 1 (t,v)-Q 2 (t,v))dv . 

The condition (15) ensures the nonpositivity of ^Q 1 (t, u). Moreover, A 2 ^ K 1 
readily implies Q 2 Q 1 . Using the induction hypothesis, we deduce that S 2 
Sj. □ 

Remark 3. Lemma 2.3 has an interpretation in terms of stochastic order or ma- 
jorization [15]. Recall that for R- valued random variables Z X ,Z 2 , the stochastic 
order ^ st is such that Z 2 ^ st Z 1 holds if and only if F(Z 2 > t) > P(Z 1 ^ i) for all 
tel. Now, let X(t) denote the Markov process representing the age of a single indi- 
vidual, defined in Remark 1, and let Y(t) denotes its number of jumps at time t, i.e., 
the "generation" of an individual. Then, Sj(t) = J2i^j Jo°° n i{ x ^)dt = P(Y(t) ^ j), 
and so, denoting by Y 1 (f),Y 2 {t) the two processes obtained in this way from 
A 1 , A 2 , we see that the properties of Lemma 2.3 are equivalent to 

A 2 ^ A 1 => Y 2 {t) > st Y^t), Vt 

A 2 < A 1 => Y 2 (t) ^stY^t), Vt . 

3. ASYMPTOTICS OF THE GROWTH RATE 

We now assume that the division rate is of the form 

K(t,x) = nip(t)B(x)l [ai+oo[ (x) 

an determine the limit of the Floquet eigenvalue of System (2) as the division rate 
tends to infinity, i.e., as the scaling parameter n tends to infinity. 
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3.1. When both ip and B are positive: proof of Theorem 2. The Floquet 
eigenvalue X F satisfies 

f &N"(t, x) + £N K (t, x) + [X F + KiP(t)B(x)l [a)+oo[ (x)] N K (t, x) = 0, 

(23) < N K (t,0) = 2niP(t) B(x)N K (t,x)dx, 

[ N K > 0, T-periodic. 

The existence of the Floquet eigenvector N K is derived in the appendix as a con- 
sequence of the Krein-Rutman theorem. In particular, the function B is required 
to be positive and of infinite integral. We normalize the Floquet eigenvector by 
requiring that 

pT poo 

(24) / / N K {t,x)dxdt = 1 . 

Jo Jo 

To establish Theorem 2, we shall think of the Floquet eigenvalue X F = Xp(tp) 
as a function of the time modulation ip. Since the division rate is a nondecrcasing 
function of ip, using the comparison principle (Theorem 1.1), we get 

X F {±) < X F (rP) < \%@) , 

where ip and ip denote the constant functions equal to the minimum or maximum 
of ip, respectively. When the function ip(t) is equal to a constant a, the Floquet 
eigenvector can be chosen to be independent of t, so that N K (t, x) = N K (x). Then, 
the Floquet eigenvalue X F (a) is actually an ordinary eigenvalue. Indeed, it is the 
only scalar (i solution of the system 

f ^N K (x)+[ f i + K aB(x)l [a , +oo[ (x)]N K (x) = 0, 

(25) < N K (0) = 2na$™ B{x)N K {x)dx, 
{ N K >0 

Integrating the latter differential equation, we get 

N K (x) = N K (0) exp {-hx-ko. f B{y)dy) . 

J a 

Using the boundary condition, and eliminating N K (0), we determine fi via the 
following implicit equation: 

p OO pX 

2 / naB(x) exp ( — fix — na / B(y)dy)dx = 1 

J a J a 

Since fi = X F (a) ^ (in fact, /i > 0, see Theorem A.l in the appendix), we have 

pQO pX 

1^2cxp(— fia) / naB(x) exp ( — na I B(y)dy)dx = 2 cxp(— fia) 

J a J a 

and, taking a := ip, we get 

X F (iP) ^ \$@) ^ ^ . 

a 

To bound the Floquet eigenvalue from below, we perform an integration by parts, 
so that: 

p oo px 

1 = 2 exp(— fia) — ft / exp f — fix — na I B(y)dy)dx . 

J a J a 
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Take now a := tp. Using the fact that < Xp(ip) and a dominated convergence 
argument, we deduce that 

/ exp ( — X F (ip)x— Kip / B(y)dy)dx ^ / cxp ( — Kip / B(y)dy)dx >0 

Ja V _ ~Ja ' Ja V ~Ja ' K ^+°° 



and so 



2exp(-A£(V0a) 



-> 1, 



K — ^-f- OO 

which shows that lim inf K-yoo X F {ip) ^ lim K _ J . 00 X F (ip) = log2/a, concluding the 
proof of the theorem. □ 
We give a numerical illustration of this theorem in Figure 5. These numerical 
simulations were produced along the lines of [3] using a monotone finite difference 
scheme. The Floquet eigenvalue was computed by applying the power algorithm to 
the one day discretized evolution operator. 



K = 10 

— k = 30 

— k = 50 

limit log 2/ a 




0.2 0.4 



0.8 1 1.2 

minimal division age a 



Figure 5. Convergence of the Floquet eigenvalue to log2/a. 



3.2. When ip can vanish: proof of Theorem 1.2. When ip can vanish, the 
conclusion of Theorem 2 fails; its proof uses the assumption that min(?/>) > 0. 
Then, we next show that more complex behaviors appear. In particular, the limit 
of X F can be a discontinuous, staircase like, function of the majority age a, and we 
shall see that the limits n — > +oo and ip — > (locally) may not commute. 

3.2.1. Intuitive derivation of the formula via a deterministic jump Markov process. 
We consider the original PDE (2), taking T = 1, 

(26) ^ 

{ n(t,0) = 2nip(t) f™n(t,x)dx. 

where ip is a 1-periodic square wave, such that if lt\ is the integer part of t, 

(27) 3 < r < 1, ip = t^ %,, r[ (i - L*J )• 

We next derive the limit of the growth rate A when k — > +oo (within this section, 
A denotes the Floquet eigenvalue A^, thought of as a function of a and ip). This 
derivation will be purely formal for the moment, leading to an ansatz the validity 
of which will be proved in the next sections. 

Intuitively, whenever the cell is old enough in the cycle (x > a) and the time is 
favorable (0 ^ t — |^J ^ r ), a division may occur. The idea is to consider that at 
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time t = 0, we have only cells of age (n(0, •) = So), and to look formally for an 
asymptotic solution 

n{t,-)~2 m ^8 a(t) , k^oo 

of the PDE (26), meaning that the population at time t consists mostly of 2 m ^ 
individuals all of age a(t), as k — > oo. If for some integer time K, we have a(K) = 0, 
so that n(K, ■) ~ 2 m ^ K ^So, using the time-periodicity and linearity of (26) we may 
expect that n(pK,-) ~ 2P m ^S , for all inte gers p ^ 1. Commuting the limits as 
k — Y oo and p — >• oo, wc arrive at the following expression for the limit of the growth 
rate as k — > oo, 

~,, , m(K) m(t) 
(28) A°°(a) := = Jim 



if log 2 t->oo<log2 

The idea is now that the trajectory (a(t))t^o, representing the age of the population, 
is produced by a degenerate (deterministic) time-dependent cadlag jump process, 
in which the state space is M + , and a particle moves to the right with unit speed, 
untill it reaches an age and a time at which division is permitted, meaning that 
a(t) ^ a and t G [0,r[+N. At every time satisfying this condition, it jumps to 
point (age is reset to zero and a division occurs). The number m(t) counts the 
number of jumps up to time t and represents the number of divisions of the cell 
which occured up to that time. 

It suffices for our argument to consider the case in which a(0) = 0. Then, the 
division times can be readily computed: if we start from age at time 0, then, 
a division should occur at every instant a, 2a, . . . ,ka as long as it is permitted, 
meaning that a, 2a, . . . , ka £ [0, r[+N. Let 

(29) K a := inf{fc E N; ka e [t, 1[+N} , 

so that K a a is the first time at which a division is not permitted, or K a = oo. When, 
for the first time, we reach a possible division time K a a that is not permitted, then 
the corresponding division occurs later, namely at the next period : [if a a] , where 
[•] denotes the upper integer part. Then, at time |~-Kad] the population divides 
(for the -Kath time) and we have n(\K a a], ■) = 2 Ka 5o- Therefore, (28) yields 

(30) e \°°{a)\K a d\ = 2 k % A oo (a) = _^_ log2 . 

I A a a 

When K a = oo, we shall adopt the convention that A°°(a) = (log2)/a. 
This is illustrated in Figure 3.2.1. 

Working with a semi-open interval [t, 1[ in (29) is consistent with the cadlag 
nature of the process a(t), however, it leads to the possibility that K a = oo, which 
is somehow unpleasant. Hence, we will use in the sequel an equivalent formulation 
of A°° , obtained by replacing [t, 1 [ by the closed interval [r, 1] , 

N a := min{fc 6 N; ka G W T } , W T := [r, 1] + N . 
Lemma 3.1. We have N a < oo. Moreover, 

(31) A°»^-^-log2 . 

\N a a\ 

Proof. If a is irrational, the finitcness of N a follows from the density of the sequence 
ka modulo 1, in [0, 1[; whereas if a is rational, we have ka G 1 + N for some integer 
k, so that N a ^ k is finite. It remains to show that the expressions (30) and (31) 
coincide. Assume first that N a a G [r, 1[+N. Then, by definition of N a , ka G]0, t[+N 
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i 88c 

division possible ; division impossible 

» 48 | 4gT 48^j 



}28b 28, ; 
1 8 a ! 

II -£ 1 time 

Figure 6. Example of the derivation of the formula. Here 2a < 
t < 3a < 1 and we obtain n(l, •) = 8n(0, •). Hence we expect 
A°°(a) = 31og(2) 

for 1 ^ fc < N a , and so, ka € [0, r[ for ^ fc < iV a , which implies that K a = N a . 
Assume now that N a a € 1+N. Then, we deduce from ka €]0, r[+N, for 1 ^ k < N a 
that (N a + k)a g]0,t[+N, for all 1 < k < N a , and so K a = oo, which by our 
convention regarding (30) leads to (31). □ 

In the remaining subsections, we show that lim K _ i . 00 A K (a) does coincide with 
the function A°°(a) (Theorem 1.2). Throughout the proof, the parameter t will be 
fixed. The proof of the theorem relies on the following strategy: 

(i) prove that the ansatz A°°(a) (staircase formula) is an upper bound for 
A K (a), for every k > 0, (Section 3.3); 

(ii) derive technical properties of the limit and of A°°(a); in particular mono- 
tonicity and right continuity, (Section 3.2.2); 

(iii) assume that a is not the right end of a step of A°° , 

(iv) consider the generational formulation of the division equation with initial 
data N K (0,x) and show that at time |~iV a a], almost every individual has 
reached generation N a and deduce the convergence for such an a. 

(v) conclude that the limit of A K as k — > oo is given by the staircase formula of 
A 00 (a), for all a > 0, using the fact that both A°° and lim^oo A K are right 
continuous functions of the parameter a. 

3.2.2. Preliminaries: technical properties of A°° and of lim K _ ! . 0O A K . We shall need 
a number of technical observations regarding the function A°° defined in (30). 

Lemma 3.2. The function a i— > A°°(a) is nonincreasing. 
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Proof. Let (ct(t))t^o denote the trajectory of the jump process defined above, start- 
ing from a(0) = 0, and let m(t) denote the number of jumps (divisions) up to time 
t (included). We have 

m(t) = sup {k G N; 30 ^ ai ^ . . . < a k < t, a i+1 - ai ^ a, a, G [0, t[+N} 

(number of divisions since time 0) 

Since the constraint a.; + i — a.; ^ a gets stronger as a increases, it is then obvious 
that m is a nonincrcasing function of a. We have, by definition of A°°, 

(32) Hm ^ = ^ = ^. 

v 7 t^+oo i [jV Q a] log 2 

Since for every i, m(i) is a nonincreasing function of a, the same is true for 

]ha t m(t)/t = A 00 / log 2. □ 

Lemma 3.3. The limit £(a) = lin^^oo \ K (a) satisfies 

£(a) = supA K (a) , 

and the map a i— > £(a), is nonincreasing and right continuous. 

Proof. Since X K (a) ^ log 2/a, and since, by the comparison principle (Theorem 1.1), 
A K (a) is a nondecreasing function of k, 

£(a) = supA K (a) < ^ . 

Furthermore, since for any fixed k, X K (a) is nonincreasing with a, so is £(a). Finally, 
since for any fixed n > 0, A K (a) is continuous with respect to a, £(a), which is a 
supremum of a family of lower semi-continuous functions, is lower semi-continuous. 
Combined with monotonicity, this leads to the convenient property: 

Va>0, £(a) = liminf£(x) = lim t(a). 

x— ¥a x— 

That is, 1(a) = lim K _ ! . 00 \ K (a) is right continuous. □ 
Lemma 3.4. The function a i— > A°°(a) is rig/it continuous. 
Proof. By Lemma 3.1, 

AU = min W T n Na =: /(a), A°» = log 2 = -^L log 2, 

|A^ a a| a|/(a)| 

where we just denote f(a) = N a a. To show that the function A°° is right continuous, 
it suffices to prove that 

< \ M 

3W = r77 m ' 

is right continuous. If g(a) ^ 1, this means that f(a) ^ [/(a)], and therefore, for 
e > small enough, we can have n(a + e) G [[nojj L na J + r [ f° r a H n < -^a an d 
N a = N a+e , f(a + s) = f(a) + N a e and \f(a)] = \f(a + e)] . Thereby, for e > 
small enough, we have g(a + e) = g(a) + rj7^yi £, which leads to the right continuity 
of g at a. If g(a) = 1, then, we can check that 



and so 



lim f(a + e) = 

S-X) + 



lim g(a + e) = 1. 

£->0 + 
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□ 

These observations leads to the following useful property (stating basically that 
A 00 is piccewisc constant) 

Lemma 3.5. For all a > 0, I a = {A°°(a;) = A°°(a)} is an interval with nonempty 
interior. Moreover, inf I a g I a ■ Finally, if we denote 

E T : {a > 0, a — sup/ a } , 

we have 

• if a g E T then X°°is continuous at a, 

• if a g E T then N a a is an integer. 

We shall see that the set E T contains exactly the right ends of steps of A°° where 
no jump occur. An infinite number of steps accumulate at the right of every point 
of E T . On Figure 1, we can see for instance that 1/2 belongs to E 2 /3 

Proof. Step 1: I a is an interval 

Firstly, this set is necessary an interval by monotonicity : if a' £ I a , then for any 
convex combination of a, a', we have 

A°°(a) = A°°(max(a,a')) \°°(6a + (1 - 6)a') ^ A°°(max(a, a')) = A°°(a). 

To prove that it is nontrivial, we use the definition of N a . By definition a, 2a, . . . (N a — 
l)a ^ W T . As W T is a closed set, its complementary is open and thus for e > 
small enough, 

}a-e,a+e[nW T =]2(a-e), 2(a+e)[C)W T =}(N a -l)(a-s),(N a -l)(a-e)[nW T = 0. 

Finally, for e > small enough, we have necessarily }N a (a — £),N a a] C W T or 
[N a a, N a (a + e)[c W T (possibly both). Thereby, for any a, there exists e > 
such that [a, a + e[ or ]a — e, a] is contained in I a which has then a nonempty 
interior (contains an interval of size e). The fact that inf I a belongs to I a is just a 
consequence of right continuity. 

Step 2: continuity on E T 

The continuity follows from right continuity of A°° and the fact that in this case, 
A°°(sup/ a ) = lirn r _j. SU p j a _o X°°(x), that is A°° is also left continuous. 

Step 3: N a a = \N a a] , forall a g E T 

Assume that a satisfies [A a a] — 1 + r ^ N a a < \N a a] . Then, for some e > 0, we 
have [iV a a] — 1 + r N a a' < \N a a~\ for any a' g]a, a + e[. Thereby, we can claim 
by definition of N a that 

N a ,^N a , \N a ,a'] < \N a a}. 

Suppose now that in addition that a g E T . By definition, as a = sup/ a , we have 
A°°(a') ^ A°°(a) for a' g]a, a + e[. Furthermore for a' g]a, a+e[, A°° takes its values 
in 

' N 

— log 2, N^N a , P ^\N a a] 
P 

therefore, for a' g]a, a + e[, we have 

and this contradicts the continuity of A°° at a. □ 
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The set E T has a last interesting property. 
Lemma 3.6. The set E T is discrete. 

Proof. Let a G E T . Thanks to Lemma 3.5, we know that A°° is constant on I a and 
therefore a' E T for a 1 G] inf I a , a[. It remains to prove the existence of e > such 
that }a,a + e[CiE T = 0. 

To prove this, let m G N \ {0} satisfy r ^ 1 — — . Choose e = 2jy 1 — . Let 
a' = a + e' G]a, a + e[. We have < e' < e. For every integer k, we have 

kN a a' = kN a a + kN a e' . 

As e' ^ 2]y~i there exists an integer k such that kN a e' ^ 1 < (k + l)N a e' . We have 
then, 

r sC 1 - - < 1 - 2N a e' < (k - l)N a e' < kN a e' ^ 1. 
m 

This leads to, adding (k — l)N a a, 

(k - l)N a a + t <(k- l)N a a < (k - l)N a a + 1. 

Since N a a £ N, this means that (k — l)N a a' € int(W T ). If we assume that a' € E T , 
we have N a ia' G N and then by construction 

{0,a',...,N a ,a'}nW T = {0,iV a ,a'} 

It follows then that 

a'NnW T = {N a ,a')N, 

which in particular implies a'N D int(W T ) = 0. This contradicts the existence of p. 
We have proved that 

] inf I a ,a + e[(~]E T = {a}, 
which ends the proof of Lemma 3.6 □ 

3.3. Upper bound on X K . We next show that A°° is an upper bound of A K , for 
each k > 0. This will turn out to be a relatively simple consequence of the heuristic 
argument, involving the earliest possible times for division, which led to the formula 
of A°°. The proof relies the generational reformulation of the division equation. Let 
a > 0, r g]0, 1[, k > be fixed. We take the eigenvector as an initial data. 

d t n t + d x 7ii + Kip(t)t[ at00 [(x)rii(t,x) = 0, 

n i+ i{t, 0) = Kip(t) J a °° ni(t, x)dx, 

n (x) = N K (0,x). 

We note by I max (t) the maximal integer i (generation) such that rii(t, •) ^ 0. At 
time t = 0, we have I max (0) = 0. Since cells need to have an age bigger than a to 
change generation, we have necessarily I max (t) ^ 1 for t < a and more generally, 

/max(<) ^ n, if t < na. 

Im*x(t)^N a if t<(N a )a. 

We introduce the notation p a = \N a a \ — 1. The integers N a ,p a are then the minimal 
integers such that 

Pa +T < N a a ^ p a + 1. 
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We use now the properties of the equation: since no change of generation may occur 
between times p a + t and p a + l,we have 

/max W < N a if p a + T < t < p a + 1. 

Now, from the definition of N K and A := A K , we have 



^2 i n i (t,x) = e A *7V K (M). 
In particular, for t < p a + 1, we have 

imax(t) iv„ 

^ 2 i n l {t,x) = ^2 i n l {t,x) = e xt N K {t, x). 

i=0 

We integrate with respect to x and obtain 

/>oo N a N a /i oo 

(33) e At / N K (t 7 x)dx = / m(t,x)dx ^ 2 Na J2 m(t,x)dx. 

JO Q JO Q JO 

But we also have from the properties of the generational formulation and the defi- 
nition Of /max, 
(34) 

Na n OO °° />00 />OC />0O 

/ riiit, x)dx = / n i (t 1 x)dx= I no(0,x)dx= / iV^O, x)cix. 
d Jo Jo jo Jo 

Combining (33-34), we get 



e At / N K (t,x)dx < 2 iv " / 7V K (0,x)cfe, i < p a + 1. 
Jo Jo 

/>oo 

As ft is fixed, we can use the continuity of / N K (t, x)dx. We have therefore 

Jo 

/>oo poo 

cxp(A K (p a + l)) / N K (p a + l,x)dx ^ 2 N " / iV K (0,a;)dx. 
Jo J0 

As N K (p a + 1, x) = N K (0, a), we have 

cxp (A K (p a + 1)) <2 JV ». 

And therefore, as expected 

A» < J-?_log2 = A°>). 

Pa + 1 

3.4. Bounding A from below. 



□ 



21 




Figure 7. Illustration of Lemma 3.7 and of Lemma 3.15: the 
trajectory of the jump process a(t) starting from a(— 1 + r) = is 
shown in bold for a = 0.22 and r = 0.6 (division is not permitted 
in the light ray regions). Here, N a = 3 and p a = 0, so X°°(a) — 
3 log 2. The minimal division age a can be increased up to a r = 
0.3, leaving A°°(a) unchanged. The conditions of Lemma 3.7 can 
be satisfied with a positive e, because, for all a' G [a, a r [, the 
combinatorial type of the trajectory of the jump process a(t) is 
unchanged (trajectory shown by a thin line, for a' close to a r ); 
choosing 8\ := 0, the times 82 and #3 can be taken to be the 
abscissae of the left of the dark grey regions of width e < (a r — 
a)/2 = 0.04. 



Step 1: Deriving inequality (12). (with e instead of r.) We derive the bound 
through the following technical lemma, in which the sequence 8\ , . . . , 9k represents 
possible successive division times of an individual of age 1 — r at time 0. 

Lemma 3.7. Suppose there exist 0i,...0/v a an d £ ^ 0, satisfying the following 
properties: 

!6>i > max(a - 1 + r,0), 
Vi<JV«-l, ff <+ i-(ff 1 +e)>a J 
Vie{l,...,iV } i+e < [^j+r, 

then the following inequality holds: 

(36) e v»(Pa+i) ^ ( 2 - e -« £ )^. 

Intuitively, this lemma shows that if there is a "tube" of positive width e, limited 
from below by such a sequence 81, . . . ,0N a and included in the set of times at 
which division is permitted, then, "enough" individuals can follow this tube, so the 
Floquct eigenvalue A K (a) can be bounded below by (iV log2)/(p a + 1), up to a 
correction of order 0(exp(— ne)) as k — > 00. 

The existence of the sequence 9\, . . . ,8jsr a will be established in Step 2 below 
when e = 0. We shall see in Step 3 that such a sequence does exist, for some e > 0, 
provided that a does not coincide with the right end a r of a step of the staircase 
function X K . This is illustrated in Figure 7. 

Proof. As for the upper bound, we consider the generational reformulation 16, in 
which K(t, x) — Ki/)(t)l[a jOC [(x) and the eigenvector N K (t — 0, x) is the initial data. 
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From the definition, we have 

5^2 < n i ( Po + l l x) = 7V K (0,a;)cxp(A K (a)(p a + l)) . 

i 

Especially, 

poo poo 

/ n l ( Pa + l,x)dx = exp(\ K (a)(pa + l)) / N K {0,x)dx. 
j Jo Jo 

We recall that N K (0,x) — for x G [0, 1 — r[. As wc look for a bound from 
below, thanks to the comparison principle, it is sufficient to consider the individuals 
that are the youngest at the beginning, that is the individuals starting at age 
(1 — r). As K{t,x) is nondecrcasing with respect to x, we know that if wc choose 
K'(t,x) ^ K(t,x) , and n! i with obvious notations, then we have thanks to the 
comparison principle 

f-OO f-OO 

V2 ! / n' t ( Pa + l,x)dx ^cxp(\ K {a){ Pa + l)) / N K {0,x)dx. 
i Jo Jo 

We choose 

!K if £ £]8 1 ,9 1 +s], x ^ a 
k if t €]6i + (i — l)e, 8i + ie], x ^ a, 
otherwise. 

From the definition of (0\, . . .0pf a ) and e, we have K'(t,x) ^ K(t,x). The most 
important property of the transition rate K' is the following. 

Lemma 3.8. Suppose that n[ are defined by the generational dynamics with tran- 
sition rate K' , then for any 1 ^ i ^ N a , 

suppn-(#i, ■) C [a, +oo[. 

Proof. We have for x < a, 9i~ x Uj]^i + (* — l) e i an d thereby K'(9i~x, ■) = 

0. Therefore for any j ' ^ 

/>oo 

n 'j+i(9i -x,0)= / A w (6>i - a;, y)n^(^ - x, y)dy = 0, 
Jo 

and using the characteristics 

n j+1 (6i,x) = n j+1 (6i - x,0) e -K K '<- e '-*+'>') da = . 

Finally, as supprip C [1 — r, +oo[ and as suppno(i, •) = suppn§ + t, we have since 
9i ^ a — 1 + r for any i 

suppn (6>i, •) C [a, +oo[. 

□ 

The latter inclusion means that for any j > 0, at time 9j, every individual is 
mature enough to divide. 

Lemma 3.9. The population ni satisfies the following equality : defining Ii{t) = 
rii(t,x)dx, we have 



Vj > 0,Vz, 



h{e j + e) = e-™I {0 j ) i 

Ii+i(8j + £ ) = e-^I l+ i{9 3 ) + (1 - e-^)h{9 3 ). 
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Proof. It is essentially based on the previous lemma. By construction, we have 
e < a. Since we already proved that supp n,i(9j, •) C [a, +oo[, we can claim 

supp rii(8j +£, •) C [0, e] U [a + e, +oo[ for i > 1, suppno(0j + e) C [a + e, +oo[ 

From the definition of the rij, we have, 

Vi,Vx ^ a+e, ti j (^+£,.t) = n^^x-^e' ^ K'{e j+ s,x-e+s)ds = n .(g j . |a! _ e ) e -«=. 
Therefore, 



ni{6 j -\- e , x)dx ~ e K£ m(9j,x)dx = e Ke Ii(9j). 
We have also, for all i ^ 0, 

n i+ i(^-+e,x)da; = f n i+ i(0,- + e - x, 0)e~ f » K'(e i+s -x+s,s)ds dx 

rij+i (9j + e — x, 0)dx 

rii+iiflj + s, 0)dx. 
Since suppnj(#j, ■) C [a, oo[ and e < a, we also have 

/•OO /"OO 

n i+ i(0j+s,O) = / K'(6j + s,x)n i (0 j +s,x)dx = / Km(9j, x)e~ RS = Ii(9j)Ke KS . 

JO J a 

We integrate and obtain finally 

n l+1 {9, + £,x)dx = (1 - e- Ke )Ii{9j). 

This completes the proof of the lemma. □ 
To conclude the proof of Lemma 3.7, we shall need the following binomial type 
representation. 

Lemma 3.10. Denote by (™J the binomial coefficients (with value if p > n), then 
we have 

Vi, h{9 3 +e) = J o (0) Q (1 - e^y(e- K y-\ 

Proof. This is readily obtained by induction from Lemma 3.9, using the fact that 
It(O j + e) = Ii(0 j+1 ), for allj. □ 
Wc finally estimate 

V2 ! / nj(p„ + l,x)dx = Y^2 l h(p a + 1). 
J o i 

It follows from the construction of K' that ii(p + 1) = Ii{9N a + £ )- Therefore, 
2 1 «J(P« + 1, x)dx = 7o(0) 2 2* (*°) (1 - OVT' - * 

= (2(1 - e- Ke ) + e - KS ) N ° = (2 - e^ 6 )*" . 
This ends the proof of Lemma 3.7. □ 
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Step 2: Construction of the sequence 9i, with e = 0. We need now to prove the 
existence of the sequence 9i used in Lemma 3.7. First we prove the following 

Lemma 3.11. There exists a sequence 9i as in Lemma 3.7 with e = 0. 

The following tj. 's define the dynamics of successive division times of a cell ini- 
tially of age x changing generation/ dividing whenever it reaches age a at a time 
which is allowed. These will be used afterwards through the construction of 

{0 if x ^ a, 
a-x if a-x [t, 1[+N, 
\a — x~\ otherwise 

\ \U(x) + a] otherwise 

We now describe the behaviour of the tk's. First of all, from the definition of N a ,p a , 
we have the following 

Lemma 3.12. The sequence ti satisfies the following properties : 

• fori< N a , tj(0) = ia, t Na (0) = p a + 1, 

• for x > a, for any i > 0, U + i{x) = ti(0), 

• x i— > ij(a;) is nonincreasing for any i, 

• Vi, tj(a:) «S ii(0) ti+i(x). 

Proof. The first point is an immediate consequence of the construction of N a ,p a . 
The second point follows from the fact that if x ^ a, then t\(x) = 0, and of the 
definition of N a . The last two points are straightforward. □ 
An important property of the sequence U is given in the 

Lemma 3.13. Suppose x > 1 — r or (x = 1 — r and A^ a a 7^ p a + 1,), i/ien we /iaue 

*JV a (a;) < Pa + 1. 

Proof. Firstly, this is obvious if x a since in this case tN a {x) = tw a -i(0). If 
.t < a, we use the following 

Lemma 3.14. Suppose x < a, let the U(x) be defined by (37). Define 

io '■= inf i. 

ti(x)£N 

Then we have 

Vn < N a , U 0+n (x) = U {x) + na, 
Vi < io, U(x) = ia — x, 
U a (x) = \i a - x] 

Proof. This follows immediately from the definition. □ 
As a consequence of this lemma, we know that if io < N a (taking n = N a — io 
in the lemma) or if io > N a , we have 

tN«(z) ^Nandtjv a (x) ^t Na (0) G N, 

therefore, in this case t]y a (x) < tjy a (0) = p a + 1. 
If io = N a , then, since x > 1 — r, 

W Q a - x<p a + l- a;<p +r. 

Therefore, we have necessarily i jv (2;) < Pa + t < p a + 1 which ends the proof of 
Lemma 3.13 and thereby of Lemma 3.11. □ 
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Step 3: Construction of the sequence 9i, with e > 0. Here, we use the notation 
a r = sup{a',A°°(a') = A°°(a)} and a/ = inf{a', A°°(a') = A°°(a)}. 

Lemma 3.15. Suppose there exists a sequence 9\, . . . , 9^ a and a r~ a ^ e ^ 0, such 
that 

6i ^ max(a — 1 + r, 0), 
Vi, C [10<J,|0<J+t], 

Vi < JV„, fli+i ^6,+e + a . 

Then there exists a sequence 8[, . . . , 9' N 

6[ > max(a - 1 + r,0), 
Vi, [fiJ.flJ + elctLflJJ.LflLI+r], 
Vi < N a , 8' 1+1 ^e' i +e' + a'. 

with e' = min(e + ^^-,r). 

Proof. We start the proof with an obvious but necessary remark : if A°° is constant 
on ]ai,a r [, then 

• e < t 

• a r — ai ^ t 

We build the sequence following this procedure : 

^=max(^j,^-^)- 

With this construction, we have 6[ ^ a' — 1 + r and since e' — 2=2. ^ e 5 

fl{ + e' < max( [^J + e', 0< + e) sC [6>iJ + r. 

Note that we also have 6[ + e' ^ 0j + e + a Z a . This helps us to check the last 
property: 

e' i+1 - e[ - a' - e' > e l+1 - a —^- -e^a'-e- a —^- = e l+1 -^-00. 

This ends the proof of the lemma. □ 
To complete the proof of the theorem, we combine the previous arguments. If A°° 
is constant on [ai, a r [, then for a < a r we can build a sequence 9i as in Lemma 3.7 
for any e < (a r — a/2) (we first build a sequence for a r — and e = and then 
use Lemma 3.15 to build a sequence with e = a r — a — 0). Using Lemma 3.7, we 
get the inequality (36) for any e < (a r — a/2). By a standard continuity argument, 
the same is true for e = (a r — a/2) = r(a, r) as defined in (12). Finally, taking the 
logarithm of (36) with e = r(a, r), we get 

, log(2-e-^) 
A (a) > A (a) — . 

This leads to (12) thanks to the following inequality, which is a consequence of the 
concavity of the function log: 

Vze [0,1], log(2-z) > (log2)(l-z) . 

This ends the proof of the theorem. □ 
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3.5. Numerical illustrations. We next illustrate Theorem 1.2 by numerical ex- 
periments, that we produced using again monotone finite difference schemes (sec 
the discussion concerning Figure 8 above). It should be noted that the bound of 
the rate of convergence given in this theorem vanishes when N a a is a multiple of 
the period T. This is confirmed by the numerical experiments: the correspond- 
ing values of a are discontinuity points of the staircase function A°°(a), and the 
convergence is seen to be slow at the left of each of these points. 



— K=10 

— K=30 

— K=50 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

minimal age for division a 



Figure 8. Convergence of the Floquet eigenvalue, when r = 1/2. 
The curves a i-> X(a) are shown for k = 10, 30, 50, 100. The limit 
A°°(a) is a staircase function. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
mininiiil age for division :i 



Figure 9. Convergence when r = 1/3. Same conventions as in Figure 8. 

Appendix A. Appendix: Existence theory for the Floquet eigenvalue 
We shall use the following notation: 

• T/j-B = x i->- B(x + h), extended to take the value zero on [0, h] if h < 0, 

• I/^ r (0,T, X) is the space of bounded T-periodic functions taking values 
in X, similarly, C por (0,T, X) is the space of T-periodic continuous (with 
respect to the time variable) functions taking values in X. 
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Theorem A.l. Assume that a > 0, k > 0, that -0 is nonnegative, bounded, not 
identically zero and that B positive, bounded, satisfying 

poo />oo 

(38) Vi, / ip(t - x)B(x)dx = / ip(t + x)B(x)dx = +oo. 

Jo Jo 

Then, there exists a unique Xf > such that there exists (iV, </>) in L^ r (0, T, L 1 (]R + )) x 
C P er(fi, T, L°°(M + )) satisfying (3-4-5). Furthermore, <p is unique. 



Before starting the proof, we give a few remarks on the hypotheses: 

• we do not need regularity assumptions on the function ip, 

i>00 

• when mm(ip) > 0, the condition (38) is equivalent to / B = +oo, in 

Jo 

the case min(?/>) = it is for instance satisfied if min(B) > or at least 
liminf +00 _B > 0, it is not optimal but other conditions would need assump- 
tions on a and k. It could be replaced for instance by 

/>oo p OO 

Vi, / Ktp(t - x)B(x)dx, Kip(t + x)B(x)dx > log 2, 



but, as we are studying asymptotic properties, we rather restrict the study 
to a case where existence does not depend on a nor on k, 

The proof of Theorem A.l is based on the method of characteristics and on the 
Krein-Rutman Theorem, as in [17], but we need more precisions in order to relax 
the regularity assumptions on ip. For the sake of simplicity we take n = 1. We 
set P{t) := J°° B(x)N(t, x)dx. Using the methods of characteristics, we have the 
following integral equations 

poo px 

N(t,0) = 2/ ip(t)B(x)N(t-x,Q)exp(-\x- ip(t - x + s)B(s)ds)dx 

J a J a 

= : CKN(;0))(t), 

poo PX 

P(t) = 2 ip(t- x)B(x)P(t - x) exp ( - Xx - / ip(t - x + s)B(s)ds)dx 

J a J a 

<t)(t,0) = 2 / ijj(t + x)B(x)cj)(t + x, 0) exp ( - Xx - / ip(t + s)B(s)ds)dx 

J a J a 

=: 4(0(.,O))(t). 

We have defined three linear operators on L^ r (0, T). These are well defined as soon 
as A > 0. Moreover, we can see L\ and L\ as operators on the space C per (0,T) of 
T-periodic continuous functions. We have 

Lemma A. 2. Under the assumptions of Theorem A.l, for any X > 0, L\ and £3 

are nonnegative, compact linear operators on C per (0,T) . 

Proof. The non negativity and linearity are obvious. We next show the compact- 
ness. If we fix A, a > then for / continuous and T-pcriodic with ||/|| ^ 1, if we 
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write g = C x (f), we have 

/>oo t*X 

g(t + h) = 2 / tp(t + x)B(x — h)f(t + x) exp (X(h-x)- ip(t + s)B(s)ds)dx 

J a-\-h J a-\-h 

o pX 

ip(t + x)B(x — h)f(t + x) exp (X(h-x)- ijj(t + s)B(s)ds)dx 



-h 



ip(t + x)B(x — h)f(t + x) exp (X(h-x)- if)(t+ s)B(s)ds)dx 

Ja+h 



- 2 

This leads to 

g(t + h)—g(t)= —2 / ip(t + x)B(x — h)f(t + x) exp (A(/i -x)-j ip(t + s)B{s - h)ds)dx 

J a J a-\-h 

poo pX 

+ 2(e Xh -l) ip(t + x)B(x-h)f(t + x)exp(-Xx- ij)(t + s)B(s-h)ds)dx 

J a J a-\-h 

poo pX 

+ 2 ip(t + x)(B(x-h)-B(x))f(t + x)exp(-Xx- ip(t + s)B(s - h)ds)dx 

J a J a-\-h 

4>(t + x)B(x)f(t + x)e- Xx [e~ £w. _ e - f; *(t+.)B(.)d,^ dXj 

= h+h + h + h. 
We deal with each Ii separately, 

(39) ^ 2||'0||||-B||e A ' l /i, 

(40) \h\^2\\n\\B\\{e Xh -l)^, 

to make J3 small, we make the following remark: for any R > a, 

poo 

(41) |I 3 | < 2\\i/>\\ \\T h B-B\\ Ll([atR]) + 4||^|| ||B|| / e~ Xx dx . 

JR 

Observe that for all a and R, \im h ^ + \\rhB — B\\ L i^ a ^ = (indeed, this is obvious 
if B is continuous, and, by density of continuous functions in L 1 ([a,i?]), the same 
is true under the present assumptions). Since A is positive, the second term in (41) 
can be made arbitrarily small by choosing R large enough, and then, the hrst term 
can be made arbitrarily small by choosing h small enough. It follows that for all 
e > 0, we have 

|/3 1 ^ e for h small enough. 
The same method can be applied to I4. Finally we have the equicontinuity of C X (B) 
where B is the unit ball for the supremum norm. Thanks to Arzela-Ascoli theorem, 
£3 is compact. 

We are now in position to apply Krein-Rutman theorem, there exist U 2 , U% nonneg- 
ative eigenvectors of C x , £3 associated to their respective spectral radii /32(A), /03(A). 
We have now 

Lemma A. 3. 

(42) PiW,Ui >0, 

(43) £ x (i>U 2 ) = p 2 (X)i>U 2 , 

(44) P3 = P2- 
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Proof. The second statement is a straightforward computation, we prove the first 
by contradiction: if U 2 vanishes, then, for some t, 

f-X 

ip(t - x)B(x)U 2 (t - 2;) exp ( - Xx - / tp(t - x + s)B(s)ds)dx = 0, 

J a 

as B > 0, it would mean thanks to T-periodicity ipB = which would lead to 
U2 = 0. Finally, as U2 does not vanish then £3(^2) > and therefore, /02(A) > 0. 
The equality comes from the duality of operators C\ and £3, we have 

P2 I i>{t)U 2 (t)U 3 {t)dt= [ £i{if>U 2 ).U 3 (t)dt = [ i>U 2 C 3 (U 3 )dt 



= P3 [ i){t)U 2 (t)U 3 (t)dt, 

JO 

therefore P2 = p 3 = P- 

To end the proof, we need to find A such that p(X) = 1. Obviously, p is a decreasing 
function that vanishes at infinity. Since 

O r-X r-X 

%p(t + x)B(x) exp {— ip{t + s)B(s)ds)dx = 2 exp ( - / i{>{t + s)B(s)ds) 

J a J a 

= 2, 

p — > 2 as A — > 0. Therefore there exists a unique A satisfying p(\) = 1. Up to a 
renormalization, (j> and P are unique, and therefore so is N. This ends the proof of 
Theorem A.l. 
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